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I. INTRODUCTION 



The meson-induced production of massive dileptons off baryons in the Drell-Yan (DY) 
process Q MB l+l-X provides a useful means of analyzing the quark structure of an 
unstable hadron, like the pion. Indeed, one can extract (see, for example, [2I, Si) in such a 
context the quark structure function of the pion and test the process independence of the 
nucleon structure function measured in deeply inelastic scattering. 

On the other hand, for large and large longitudinal momentum (carried by a pion's 
quark constituent) the hadronic differential cross section for the production of a massive 
lepton pair via the annihilation of an antiquark and a quark in the colliding hadrons — 
pion and nucleon, respectively — involves the pion distribution amplitude (DA) in order to 
describe the pion bound state. Tuning the DY reaction to the kinematic edge of the phase 
space, where the antiquark u from the pion is far off-shell, i.e., — >■ 1, it is sufficient to 
treat the quark m, originating from the nucleon, as being nearly free and on-mass-shell, so 
that the bound-state details of the nucleon become irrelevant In those circumstances 
the process MB l^l~X reduces to tx^u — »• 11^ n'X and the corresponding amplitude 
becomes calculable within perturbative QCD. The binding effects of the pion state are taken 
into account by means of the pion DA (which is the pion wave function integrated over 
transverse momenta), making this type of process suitable to test the details of proposed 
nonperturbative models for the pion DA. 

Moreover, the angular distribution of the produced lepton pair (actually the /i+), relative 
to the pion direction, depends in a sensitive way on the pion DA. Hence, measuring the 
angular distribution parameters A, /i, v (for an unpolarized target) and /i, p (for longitudinally- 
polarized protons), one can extract useful information on the shape of the pion DA 0, g. 
This is particularly important with regard to a proposed experiment by the COMPASS 
collaboration to collect high-precision data from the scattering of a pion beam off a polarized 
target 0]. 

In the present work we will consider (i) the inclusive production of dimuons from the 
hard scattering of pions on an unpolarized nuclear target and (ii) an analogous process with 
longitudinally polarized protons. The appearance of a single transverse-spin asymmetry in 
the latter case is related to an imaginary part, which, may even dominate the dimuon angular 
distribution 0. Because this contribution is very sensitive to the pion DA, the single-spin 
asymmetry can be used in conjunction with the experimental data in order to select the 
most preferable pion DA. Our attention is focused in both cases on the role played by the 
pion bound state in terms of the pion DA. 

The paper is organized as follows. In the next section, we present the theoretical back- 
ground of the DY model in vrA^ collisions, taking into account the valence-quark bound state 
of the pion, which we model with the aid of nonlocal QCD sum rules l8[ in comparison with 
the Chernyak-Zhitnitsky 0| model and the asymptotic pion DA [iriill, [13 • Section [ml 
contains the presentation of our results for the angular distribution parameters A, /i, i/, and 
comparison with the data from the E615 experiment at Fermilab [s^. In this section, we also 
make predictions for the angular parameters p, and z/, using the aforementioned pion DAs. 
Particular attention is devoted to the endpoint behavior of different pion DAs by considering 
an azimuthal moment of the pion DA, discussed before in 0. Finally, our conclusions are 
provided in Sec. HVl 
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II. DRELL YAN PROCESS OF DIMUON PRODUCTION 

This section describes the theoretical method used to calculate the angular distribution 
of the /i+ in the DY process in a nN collision in terms of the parameters A, /i, v (unpolarized 
target) and /i, p (longitudinally polarized protons). 




FIG. 1: Graphical representation of the Drell-Yan process 'it~N — > fi~^^~X. Symbols p = Pt^, 
P = Pn and q = —x^P + XuP + qr-, where Xu = x.,^, and x^ = x^ (see for details Sec. Ill Ap . mean 
four-momenta. Here, and below, solid lines denote leptons and also the quarks emerging from the 
ttN collision, double lines indicate bound states (with their corresponding wave functions being 
illustrated by shaded blobs), whereas the (red) curly line stands for the exchanged hard gluon and 
the (blue) weavy line represents the highly virtual photon. 

The DY process is the dominant mechanism to produce lepton pairs with a large invariant 
mass in hadronic collisions, like vr^A^ scattering. In the context of this model a massive 
lepton (muon) pair is created through the electromagnetic annihilation of an antiquark from 
the beam pion and a quark from the nucleon target, as depicted in Fig. [TJ Moreover, in 
the kinematic region of large a;^ = — »■ 1 (or xl — > 1, see next section), the antiquark in 
the pion is subject to bound-state effects encoded in the valence-quark pion DA, while the 
quark from the nucleon is nearly on shell. Then, the pion bound state can be resolved via 
hard-gluon excharige between the annihilating antiquark and the spectator in convolution 
with the pion DA j5||, displayed diagrammatically in Fig. [2J This provides a means to test the 
validity of proposed models for the pion DA, derived from nonperturbative QCD calculations, 
because the angular distribution of the produced muon pair depends in a sensitive way on the 
shape of the pion DA. Contributions from higher-order hard-gluon exchange are suppressed 
by powers of as and will be ignored; those due to evolution will be taken into account in 
leading order (LO). Alternatively, ignoring chromodynamic binding effects, one may extract 
in the context of the DY model, the structure functions of the pion and the nucleon [1, 0] 
or focus on quark density functions j3, These issues are outside the scope of the present 
investigation. 

To continue with our quantitative analysis, we first present some explanations on the 
kinematics of the DY process together with the definitions of the dynamical parameters 
which describe the angular distribution in the hadronic differential cross section. 
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FIG. 2: Leading (perturbative) contributions to the QCD-amplitude describing the Drell-Yan 
process 7r~A^ fi~^fj,~X at large and xl, the latter defined in Eq. (j2.9p . The designations are 
the same as in Fig. [TJ 



A. Kinematics and variables of the DY process 



For the convenience of the reader, the relevant kinematic parameters and dynamical 
variables for the DY process MB I'^l" X (with M, B denoting a meson and a baryon, 
respectively) are compiled below in conceptual groups and in conjunction with Fig. [1) 



Momenta 

• = p: Pion momentum. 

• pn = P: Nucleon momentum. 

• Pu'- Momentum of the annihilating antiquark u emerging from the pion, p^ = Xup. As 
Xu ^ 1 (i.e., u far off-shell), becomes large and far spacelike. 

• Xu'- Longitudinal momentum fraction (light-cone variable) of the annihilating u from 
the meson. In the kinematical regime, we are considering, x^ = x^^. 

• < k.^ >: Average of the square of the transverse momentum of the annihilating u from 
the meson (kf, -C Q^). 

• Pu'- Momentum of the annihilating quark u emerging from the nucleon, p^ = XuP- 
Because — > 1, it is sufficient to consider the w-quark to be nearly free and on-shell: 
Xu = xn (quark masses and transverse momenta neglected). 

• Q = m^^: Invariant mass of the massive lepton (muon) pair, or, equivalently, momen- 
tum of the virtual photon created by the annihilated quarks (Q^ = q^)- In this analysis 
we consider values in the range = (4 — 81) GeV^. 

• s: Squared invariant mass of the initial hadrons, i.e., s = {p + P)^. In our numerical 
analysis we use s = 100 — 400 GeV^, which covers the 252 GeV vr" beam of the E615 
Collaboration at Fermilab [sf. 

Masses 

• m,r: Pion mass. Because s ^ m^, we set m.^ = 0. 
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• m: Bare quark mass. Because s ^ m^, we set m = 0; 

• rrtTv: Nucleon mass. 
Kinematic Variables 

• T = Q'^/s: Scaling parameter. 

• p = Qt/Q'- This scaling parameter is a measure for the squared transverse momen- 
tum Q"^ = —q^ of the virtual photon in the hadronic center-of-mass frame (c.m.f.). 
Inspection of Fig. [1] reveals that appears only in combination with x^p, namely, 
Qt + XuP- Then, in order to separate both terms unambiguously, one has to demand 
that gx should not contain any part of the pion momentum p, cf. Eq. (12. 6p . 

m xl = 2(5^/1/5 < 1: Longitudinal momentum fraction (Feynman x) of the lepton pair 
(associated with the virtual photon) in the hadron c.m.f. (Q^ = and Qt are 

the photon-momentum components parallel and perpendicular, respectively, to the 
incident pion momentum in the hadron c.m.f. 

and Neglecting the quark transverse momentum and 

mass (as s becomes very large), the quantities and xjq can be identified with 
Xu and Xui respectively. Combining the two equations above, one finds [1] x.„^i^ = 

±xf + (4 + 4r)'/'] /2. 
Angular distribution parameters 

• 6: Polar angle measuring the fi'^ direction with respect to the t-channel (or Gottfried- 
Jackson system of axes), i.e., cos 6* = p^ ■ p^. The definitions of other choices of axes 
(frames) can be found, e.g., in 

• 0: Azimuthal angle between the massive lepton-pair plane and the plane of the incident 
hadrons in the lepton rest frame. 

• A, fi, v: Angle-independent coefficients, depending on and x/,. These three param- 
eters control the dilepton angular distribution and are sensitive to the shape of the 
pion DA. 

• /i, v: Angle-independent coefficients, depending on and x^, induced by the imagi- 
nary part to the DY amplitude in the polarized case. They are sensitive to the shape 
of the pion DA. 

In our analysis we consider values of s much larger than the nucleon mass, s > 100 m^, 
so that this be neglected, though we derive the exact result with 7^ and proceed then 
with the approximate expression with rriN = 0, where it is applicable. The momenta p and 
P are on the light-cone and, hence, we have 

p^ = 0, P^ = m% ^0, 2{p- P) = s -m% = § ^ s . (2.1) 

For the cross-section calculation, one can appeal to the optical theorem and set the upper 
right line in Fig. [1], denoting the ci-quark, on mass shell. Then, 

2{p.q^) = -^ = ^ (2.2) 

1 Xu -L Xu 
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and, due to g = (x„ + Xu)p + Qt, one obtains 



1 + 



S Xu I 1 -|- 



Xn, S 



Xu (s - m%) 
2 n2 



or 



1 + 



1 Xu 
P' 

1 -^M 



1 + 



2 ^]\f 



(2.3) 



r 1 



Xu S 



(2.4) 



which fixes the value of Xu as a function of x^, p and r. We also see that a non-zero nucleon 
mass is important only at very small Xu — mj^/s. At this value of Xu the value of Xu is close 
to 1, so that, neglecting rriN, one finds 1/ [1 + {m%/Q'^)/{l + p^)]. Hence, the difference can 
be sizable and, especially at ~ 2m% and p ~ 1, it can reach as much as 20%. 

In the following considerations (and, in particular, in the analysis of the polarized DY 
process), we will use the longitudinal-momentum fraction of the photon, x^. To this end, 
recall that in the hadron c.m.f. one has 



p=^(l,0,0,+l), P= ^(1,0,0,-1), gx = (0,gxi,g±2,0) 



(2.5) 



and, using Eq. 02.21) . one finds 

q± — aT P , At 



1 -^u 



(l± + (lL, Ql = XuP+ {xu -aT)P = — (xo,0,0,xl) 



with xl being defined by 

^0 = Xu -\- Xu — ClT } = Xu — Xu -'r ClT ■ 

It is convenient to recast xl in terms of Xu, p, and r to read 



_ l + p2 

Xu • 



X-ji 



(2.6) 
(2.7) 

(2.8) 
(2.9) 



The inverse relation 



XO: 



Xi + v^^|T4(rTp2)7 



(2.10) 



defines X rfjj as a function of x^, p, and r. 
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B. Drell Yan vrA^ process with a pion bound state 



There have been attempts [sl and, in particular j^, to test the compatibihty of different 
pion DAs with the DY angular distribution measured in the scattering of pions off protons. 
The first work suggested that the CZ pion DA (and some variant of it) fits the unpolarized 



data [3] better than the asymptotic DA, or narrower convex versions [13||. The second 
work focused on the possibility of longitudinally polarized protons as a target and discussed 
additional angular parameters that vanish for an unpolarized target. It was argued that 
future polarized experimental data on these processes may be able to distinguish among 
various pion DAs because of the high sensitivity of these parameters to the particular shape 
of the pion DA. Though there is still no such data available, meanwhile important knowledge 
has been collected that is rather unfavorable for the endpoint-dominated type of two-humped 



pion DAs, like the CZ one. For instance, new high-precision lattice simulations |14j . Il5l . Il6l | 
give a new level of detail for the second Gegenbauer coefficient 02, yielding values around 
0.2 to 0.24 at a momentum scale of 2 GeV, well within the range suggested by the nonlocal 
QCD sum-rule estimates [s^ — hereafter referred to as BMS — and supported by the analysis 



in 13, 18, 19 



of the CLEO data 2^ on the 71 — 7 transition form factor using light-cone 
QCD sum rules. This 02 value is about two times smaller than its counterpart of the CZ pion 
DA and cast serious doubts about the consistency of this model with experiment. Moreover, 
the aforementioned BMS CLEO-data analysis has confirmed the earlier Schmedding and 



Yakovlev j2]| findings which excluded the CZ pion DA at least at the 2a level. On the 
other hand, the asymptotic pion DA seems to be also excluded, given that its Gegenbauer 
coefficient 02 is identically zero and the CLEO-data analysis relegates this DA outside the 2a 



error ellipse (for the most recent rigorous analysis, see [l9j), while being also incompatible 



with the lattice results [IJ, Il5|, Il6[ . Further theoretical arguments and details can be found 
in j22|. See also [23J for a recent development of the nonlocal QCD sum- rule approach and 
the extraction of the pion DA. 

Therefore, it would seem reasonable and timely to upgrade the calculation of the DY 
angular parameters of the 7r~N hard-scattering process by taking into account the recent 
developments quoted above. To continue, we first recall the definition of the pion DA, 
ipTr{y, (Iq), which specifies the fractional longitudinal momentum y of the valence-quark con- 
stituents in the pion at the normalization scale /Xq. At the (leading) twist-two level it is 
defined by the following matrix element 



(0 I J(;2)7'^75C(z,0)m(0) I 7r(P)) 



22=0 



iP'-f^ [\ye'y^^'">V^{y,fil) , (2.11) 
Jo 



where = 130.7 ± 0.4 MeV [2J] is the pion decay constant defined by 



(0|rf(0)7;.75M(0)|7r+(P))=zp^^. 



(2.12) 



(2.13) 



Above, a straight path-ordered Fock-Schwinger connector j25| (Wilson line) C(0, z) = 
V exp[—igsj^ t"-A^(y)dy'^~\ has been inserted to preserve gauge invariance of the operator 
product. In the following, we use the light-cone gauge which reduces the contribution of the 
Wilson line to unity. The normalization scale, /ig, of the pion DA is related to the ultraviolet 
(UV) regularization of the quark-field operators on the light cone in fl2.12p . whose product 
becomes singular for z"^ = 0. 
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The pion DA cannot be derived from first principles, but has to be inferred from non- 
perturbative QCD models (or from experiment). This not withstanding, its evolution is 



governed by perturbative QCD (lol . 



121 ] and can be expressed in the form 



ds 



(2.14) 



where ipTr{s,fil) is a nonperturbative input determined at some low-energy normalization 
point /Xg ~ 1 GeV^, where the local operators in Eq. (12. lip are renormalized, while 
f/(?/, s; /i^, /ig) is the evolution operator from that scale to the observation scale /i, calcu- 
lable in QCD perturbation theory. In the asymptotic limit, the shape of the pion DA is 
completely fixed by pQCD to be (p^^^{y) = Qy{l — y) [l3, 11, 12], with the nonperturbative 
information being solely contained in the pion-decay constant f^^. 

In the leading-twist approximation of the pion DA, in which we are working, y9,r(|/,/^o) 
can be expressed in terms of the Gegenbauer polynomials, which form an orthonormal set 
of eigenf unctions. Then, 



l + a2(/x^)Cf (2y 



+ a,ifj,l)Cfi2y 



1) + 



(2.15) 



with all nonperturbative information being encapsulated in the expansion coefficients a„. 
Depending on the nonperturbative approach applied, these coefficients can be calculated via 
QCD sum rules 0, QCD sum rules with nonlocal condensates jsj, or lattice simulations 
HESS. More details can be found in the original papers already cited, while we here 
restrict ourselves to the obtained results, which we quote in Table [H The shapes of these 
DAs are displayed in the left panel of Fig. [3l whereas the evolution effect is illustrated in 
the right panel of this figure in terms of the BMS pion DA. Note that the underlying quark 
virtuality employed, is Ag = 0.4 GeV^ [13] corresponding to a correlation length of the scalar 
quark nonlocal condensate of about 0.31 fm. As one anticipates from this figure, the key 
characteristic of the BMS-type pion DAs is that the endpoint region a; — > 0, 1 is strongly 
suppressed — not only relative to the CZ pion DA, but even with respect to the asymptotic 
one (for mathematical details, see [2^). In contrast, the double-humped shape — reminiscent 
of the CZ pion DA — turns out to be of minor importance 27 . 



TABLE I: Pion DA models used in the analysis. 



7r-DAs]^2=i Gev2 


asymptotic 


BMS [8] 


BMS "bunch" [8,18\ 


CZ [9] 







0.20 


[0.13,0.25] 


0.56 







-0.14 


[-0.04,-0.22] 





higher 





negligible 


negligible 






Let us now make some remarks on the evolution of the pion DA in LO of perturbative 
QCD. Taking into account only the first two Gegenbauer coefficients, one obtains 



¥.^0(l/, f^l) = 6y{l -y)[l + 40 (/x^) (2y - 1) + (^|) (2y - 1) 
where 0.2^ (fip) and a\'~'{fJ'p) are given by 



(2.16) 



as(/i|) 
as(/io) 



7i°V{2feo) 



(2.17) 
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2.5 



1 . 5 



0.5 




'^TT {y, M^) 




0.4 y 0.6 



FIG. 3: The left panel shows the "bunch" of pion DAs (green broken lines), derived from nonlocal QCD 
sum rules (for a summarized exposition, see, e.g., 22]), in comparison with two extreme alternatives: the 
asymptotic DA 10, 111, — dotted line — and the CZ model Q — red long-dashed line — at the momentum 
scale /x^ w 1 GeV^. The green solid line inside the "bunch" represents the BMS model The right panel 
illustrates the effect of one-loop evolution on the pion DA, exemplified by the BMS model, in comparison 
with the asymptotic solution (continuous convex line). The double- humped solid line represents V'bms(-^) 
1 GeV^, while the broken lines mark (p]^'^g{x) at 4, 20, and 100 GeV^ (with the larger scale corresponding 
to the larger value of the DA at the middle point). 



The expressions for the anoinalous dimensions and the beta-function coefficient are 
hsted, for example, in [17|. One sees from the right panel of Fig. [3] that the effect of the 
inclusion of the LO diagonal part of the evolution kernel is indeed important. This figure 
shows how the BMS pion DA v^bms(^) evolves from the normalization scale of 1 GeV^ 
(double-humped solid line) to higher momentum values at 4, 20, and 100 GeV^ (broken 
lines), with the larger scales corresponding also to larger values of the DA at the middle 
point. The asymptotic profile (continuous solid line) is displayed for comparison. 



C. DY reaction with an unpolarized target 

The angular distribution of the yU"*" in the pair rest frame can be written in terms of the 
kinematic variables A, /i, u as follows 



rfV(7r- + N + ij,- + X) 

dQ'^dQ'^dxL d cos 



oc N{x, p) ^1 -|- A cos^ 9 + iism 29 cos i 



U 2/1 

H — sm 9 cos 
2 



where la 



A(5,p) = ^|(l-5)'[(Im/(£))2 + (F + ReJ(x))2] -(4 
-^Pxf|(1-5) [F + Rel{x)]+ p^xF^ , 

o 

i/(x,p) = -_p2^ (1 -x) F [F + Re/(5)] , 



2\ 2 ~2 172 
— P ) P X F 



p (x, p) 



N{x,p) 



2 I (1-5)^ [(Jml{x)f + (F + Re/(5))^] + (4 -f p^) p^ 5^ i^aj 



(2.18) 

(2.19) 
(2.20) 
(2.21) 
(2.22) 
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FIG. 4: Angular definitions of tlie Drell-Yan process in the center of mass frame of the produced 
massive lepton pair. The axis zqj denotes the pion direction in the Gottfried-Jackson (GJ) frame. 



with 



XL+ \/xl + 4:{l + p'^)t 
2(1 + P^ 

Also displayed is the normalization factor of the cross section, A^. The abbreviations 



x{xL,p) = , . (2.23) 



F= [\y^^ (2.24) 

Jo y 

Jo y [y + x-l + ie) 

axe functionals of the pion DA and, therefore, depend on the evolution momentum scale 
Q'^ Q'^. Note that the inverse moment (12.241) of the pion DA plays a crucial role in the 
description of several form factors of the pion in perturbative QCD [27|] . We will have more 
to say about the choice of the evolution scale in the next subsection. Note also that the 
denominators in Eqs. (I2.24p and (12.251) originate from the gluon and quark propagators in 
the subprocesses shown in Fig. [2], respectively. Before presenting results for these angular 
coefficients, let us first discuss what changes are induced when the nuclear target is polarized. 



D. DY reaction with longitudinally polarized protons 

In the polarized DY process the angular distribution of the /x^ contains two additional 
parameters /2 and i/, entering additively Eq. (12.181) with the same angular structure as [x and 
I/, respectively, being, however, proportional to the target longitudinal polarization S£ : 

- -27V sepxF(p{x,Q^) 

^ (1-x)' [(F + Re/(i))' + 7rV(5)'] +(4 + p2)p2x2F2''™''' 

iy{x,p) = 2p/i(£,p) . (2.28) 

where Xp = t/x, /j,'^ is the evolution scale for the nucleon parton distributions. Note that 
all momenta refer to the hadronic c.m.f. The polarized parton distributions used in our 
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analysis are taken from [2^, whereas for the unpolarized structure functions we use the 
parameterization of Ref. 29|. To evolve these distributions from their normalization scale 
/Xq = 4 GeV^ to the scale fi^ = = 16 GeV^, we employed the Fortran codes supplied by 
the authors of these papers on the Durham web site.^ 

Lacking radiative corrections to the process under study, we cannot fix the evolution scale 
unambiguously. In an effort to get a measure for the entailed uncertainty, we have analyzed 
the dependence of our results on the choice of the evolution scale by varying fi"^ in the range 
/i^ = and /x^ = 2Q^. We found that, depending on the model pion DA used, the 

variance of the calculated p, parameter lies between 2% and 17%, with moderate sensitivity 
to the adopted p value. Details are given in Table [TTl for the (large) xl interval, in which 
the one-gluon exchange is still a good approximation and only the pion bound-state is of 
importance. Recall in this context that the only evolution effect in the case of the asymptotic 
pion DA stems exclusively from the DGLAP evolution of the nucleon parton distributions. 
On the other hand, in the case of the BMS and the CZ pion DAs, the total evolution effect 
is the result of the combination of the ERBL evolution of the pion DA and the DGLAP 
evolution of the nucleon parton distributions (unpolarized and polarized) in p^^ci (cf. Eq. 



TABLE XL Change of predictions for fj,{xL,p), shown in Fig. [71 due to different settings of the 
evolution scale. 





vr-DAs 


asymptotic 


BMS [8] 


CZ [9] 


p = 


0.06 and 0.5 < xl < 0.84 


2.0 - 2.5% 


12 - 2.5 % 


17-2.5 % 


p = 


0.30 and 0.5 < xl < 0.83 


2.0 - 2.5% 


9-0 % 


13 - % 


p = 


0.50 and 0.5 < xl < 0.81 


1.5 - 2.0% 


8-0 % 


13 - 8 % 



It was noted in [61] that the angular moment of the pion DA, defined by 



M 



ang 



sin 29 sin da{s£ = 1) = —2 ir p x F ip{x,Q ) p 



nucl 



(2.29) 



is particularly sensitive to the x (or xl) endpoint region and can be used in comparison with 
experimental data in order to distinguish pion DAs which behave differently exactly in this 
region. We have, therefore, included predictions also for this quantity (see next section. Fig. 
|8]). In similar context, it is important to consider the (experimental) single-spin azimuthal 
asymmetry (SSA) 



A 



da{si = +1) — da{si = —1) 



da{si = +1) -|- da{si = —1) 
after averaging the cross sections over the polar angle G [0, vr] 

p p{se = +1) sin 20 



A{(p,XL,p) 



2 + A + I z/ cos : 



(2.30) 



(2.31) 



In Eq. (12.301) . cr(s^ = +1) and a{s£ = —1) denote opposite hehcity states of the longitudi- 
nally polarized target. Notice that in order that the longitudinally polarized nucleon parton 



^ http://durpdg.dur.ac.uk/hepdata/grv.html — for GRV95 and 
http://durpdg.dur.ac.uk/pdfiib/gchrmann/pdf/wclcome.html — for GS96. 
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FIG. 5: Results for the angular distribution parameters A, u, and as functions of = for 
different values of p = Qt/Q- Predictions for the Lam~Tung [30[] combination, 2v — {1 — X), are also 
displayed. The green strip contains the results for the pion DAs calculated with nonlocal QCD sum 
rules 
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18l |. The (blue) solid line corresponds to one of these endpoint-suppressed DAs, termed 



BMS, while the dotted (black) solid line shows the result for the asymptotic DA, and the (red) 
dashed line is the prediction for the endpoint-dominated Chernyak-Zhitnitsky DA 0]. One-loop 
evolution of the pion DAs to each maesured value is included. The data were taken from [sj 
and were evaluated as explained in the text. 



distributions can transfer its polarization to the azimuthal distribution of the massive lepton 
pair, one needs an imaginary part and the interference of amplitudes with a phase difference 
between them. Both ingredients are provided here by the pion DA and the hard-gluon ex- 
change. Ignoring pion bound-state effects in the treatment of the DY process, one would have 
to include in the hard cross sections perturbative QCD radiative corrections proportional to 
as in order to create a SSA. 
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FIG. 6: Angular parameter fl{xL,p) as a function of xl for three different pion DAs, evaluated at 
three different values of the scaling parameter p using one-loop evolution. Results are shown for 
the asymptotic pion DA [lO , 11 , 12i] (black dotted line) , the BMS model 0] (blue solid line) and the 
CZ one ^9|] (red dashed line). The upper row corresponds to a center-of-mass energy s = 100 GeV^, 
whereas the lower one has s = 400 GeV^. 



III. COMPARISON WITH EXPERIMENTAL DATA 



In this section we present our results for the angular distribution parameters A, /i, (versus 

^ Xu) for the unpolarized DY process and compare them with the available experimental 
data. We also include predictions for the parameter ft (versus xl), which is nonzero only 
in the polarized DY process, and the SSA A. 

All results were obtained using different models for the pion DA, already mentioned, and 
they are shown at different values of Qt, i.e., at different values of the scaling parameter p. 

To be specific, we compare in Fig. [5] our theoretical predictions for Xi^iV with experi- 
mental data from E615 [3| in that x^^ region reported by this collaboration using a 252 GeV 
7r~ beam interacting in an unpolarized tungsten target. 

We used for convenience the Gottfried-Jackson frame and included in our analysis the 
data sample listed in their Table VIII. Because the range of the probed transverse photon 
momentum, Qt-, is in our opinion too large for averaging, we followed another strategy 
than the authors in Ref. [sj. Notably, we adopted some value of the scaling parameter 
p = Qt/Q and required the momenta Qt and Q to be within the reported window of the 
measurement in [3|. We evaluated this way the angular parameters shown in the figures. 



obtained with different pion DAs, including one- loop ERBL evolution with = m^^^, and 
for three different values of p = 0.06, 0.3, 0.5, whereas s = 500 GeV^ in accordance with [3*]. 
Strictly speaking, the results at too low values of p < 0.06 are, actually, not compatible with 
factorization because then Qr becomes of the order of Aqcd- They are shown here merely 
for illustration purposes. 

One important observation from Fig. [5] is that the parameter /i increases proportionally 
with p, as it is qualitatively expected from Eq. (12.201) (though there is an additional term 
proportional to p^, absent in the parton DY model). This increase turns out to be small for 



the asymptotic DA [10|, lll| 
while for the BMS strip 




whereas it becomes rather too strong for the CZ model [9| 
this enhancement with p is moderate and provides best agree- 
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FIG. 7: Evolution effect on the angular parameter nixj,, p ) v s. x r in several steps of the scaling 



parameter p. Results are shown for the asymptotic pion DA 13, U, 12] (left panel), the BMS model 
DA fs'] (central panel) and the CZ one [o] (right panel). The following designations are used: The 
solid lines within each color group correspond to the smallest p value. The short-dashed lines denote 
the results for the intermediate p values, whereas the long-dashed lines represent the results with 
the largest p values. The group of the blue lines covers the range p = 0.001, 0.003, 0.006; the group 
of black lines gives the results for p = 0.01,0.03,0.06, whereas the red lines are associated with 
the values p = 0.1,0.3,0.5. The upper row corresponds to a center-of-mass energy s = 100 GeV^, 
whereas the lower one has s = 400 GeV^. 

ment with the data. A similar behavior is seen also for u in this figure, which, according 




FIG. 8: Angular moment TWang as a function of for three different pion DAs, evaluated at three 
different values of the scaling parameter p using one-loop evolution. Results are shown for the 



asymptotic pion DA [10|, lUl, llM] (black dotted line), the BMS model [g] (blue solid line) and the 



CZ one 0] (red dashed line). The upper row corresponds to a center-of-mass energy s = 100 GeV^, 
whereas the lower one has s = 400 GeV^. 
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to Eq. (12.211) . should increase proportionally to . Also here the BMS strip compares most 
favorably with the data relative to the other options. On the other hand, the Lam-Tung 
combination 2z/ — 1 + A jsO], which is the analogue of the Callan-Gross relation in deeply in- 
elastic scattering (and also reflects the kinematical nature of the azimuthal asymmetry |3l[|). 
is badly violated by the data. 

This trend is in agreement with the theoretical predictions above approximately Xu ~ 0.6, 
though, at very large close to the kinematic limit, all tested model DAs tend to fall stronger 
than the data. 

Turn now attention to the polarized case. Figure E] shows the predictions for the angular 
parameter /2 as a function of xl for three different values of p = 0.06, 0.3, 0.5 and two center- 
of-mass-energies s = 100 GeV^ (upper row) and s = 400 GeV^ in the range expected to be 
covered in the COMPASS experiment. As before, the results shown were obtained for the 
asymptotic, BMS, and CZ pion DAs. [Similar results (not shown) hold also for the parameter 
ly, which is interlinked with p.] 

The result of LO evolution of the pion DA as well as of the nucleon structure functions 
in the angular parameter p is the detailed image of curves shown in Fig. [71 

To unravel the behavior of the various pion DAs in the endpoint region, we also examine 
the angular moment (cf. Eq. (12.291) ) and display the results in Fig. [HI As mentioned in the 
previous section, this quantity was proposed in as a sensitive measure for the angular 
modulation of the pion DA at large xl- One sees from this figure that this seems indeed to 
be the case, especially for not too small p values and, once experimental data will become 
available, it could potentially serve to discriminate among different pion DAs. 

To complete the analysis, we include 3D-plots of the parameter /2 as a function of xl and p 




FIG. 9: 3D-plots of the angular parameter p{xL,p) for three different choices of the pion DA: 
Asymptotic 10, 11, 12 1 (left), BMS model |^ (center), and CZ model (right). The upper row 
corresponds to a center-of-mass energy s = 100 GeV^, wher eas the lower one has s = 400 GeV . 
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FIG. 10: 3D-plots of the azimuthal asymmetry A{(j),XL,p = 0.3). All assignments are as in Fig. [9l 



(Fig. [9]) and analogously 3D plots of the angular asymmetry A{4>, xl, p) versus xi and (Fig. 
[To!) , the latter for p = 0.3. Note that both quantities have been evaluated for two different 
values of the center-of-the-mass energy s = 100 GeV^ (upper rows) and s = 400 GeV^ (lower 
rows), as expected for the COMPASS experiment 0]- 

The observed increase of A with xl suggests that the single-spin asymmetry of the muon- 
pair angular distribution is associated with the valence Fock-state contributions in the pion 
DA. This behavior is valid for all considered pion DAs and does not significantly depend 
on the value of the scaling parameters p and r, though its size decreases with increasing s. 
We have verified that the results do not change significantly for larger p values. If future 
experimental data would confirm such a single-spin asymmetry, this would point to a new 
mechanism for generating a new nontrivial phase in QCD factorization leading to T-odd spin 
asymmetries that mimics the effect of a true T violation. 

IV. CONCLUSIONS 

Our objective in this work was (i) to update previous results on the unpolarized 
Tc^N —>■ fi^p,~X DY process and (ii) to make detailed predictions for the angular distri- 
bution parameters for the hard-scattering of pions on longitudinally polarized protons. The 
single-spin asymmetry, predicted here for various pion distribution amplitudes, may soon 
become amenable to experimental check at COMPASS. To this end, we have presented an 
updated analysis of the DY process with the inclusion of the pion's DA, the latter based on 
the theoretical appraisal of the theoretical situation obtained within the context of nonlocal 
QCD sum rules and supported by recent high-precision lattice calculations and other exper- 
imental data from the pion-photon transition. Though the existing data on the unpolarized 
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DY Tt~N process cannot single out one particular pion DA with little ambiguity, 

the "bunch" of the BMS DAs, derived from nonlocal QCD sum rules, seems to comply most 
favorably with the E615 data. Given the distinctive behavior of all considered pion DAs 
with respect to the longitudinal momentum fraction, carried by the annihilating quark from 
the pion in the polarized DY process, one may hope that measuring the angular moment 
in the planned COMPASS experiment may lend quantitative support for one or the other 
proposed pion DA. 

On the other hand, it is also important to consider another (and actually the most com- 
mon) mechanism of the tt + N Drell-Yan process, treating the pion structure in terms of 
parton distributions rather then the pion DA. This mechanism is the dominant one at moder- 
ate values oi xl- There are also related contributions to the SSA due to the imaginary phases 



emerging in the short-distance subprocesses [3^, |33|. These contributions have recently been 



studied for the kinematics of RHIC and J- PARC [3J|. The investigation of the relevance of 



this mechanism in the case of the COMPASS kinematics requires further investigation. 



Acknowledgment s 

We would like to thank A. V. Efremov and S. V. Mikhailov for discussions and useful 
remarks. We are grateful to R. Bertini, F. Bradamante, and O. Denisov for discussions 
on the COMPASS experiment and its potential extension to the Drell-Yan process. Two 
of us (A.P.B. and O.V.T.) are indebted to Prof. Klaus Goeke for the warm hospitality at 
Bochum University, where part of this work was carried out. N.G.S is grateful for support 
to BLTP@JINR, where this work was completed. O.V.T. is also indebted to M. Anselmino, 
F. Balestra, R. Bertini, A. Kotzinyan, and G. Pontecorvo for the warm hospitality and useful 
discussions at INFN (Torino). 

This work was supported in part by the Deutsche Forschungsgemeinschaft, grant 436 
RUS 113/881/0, the Heisenberg-Landau Programme (grant 2007), the Russian Foundation 
for Fundamental Research, grants No. 06-02-16215 and No. 07-02-91557, and the Russian 
Federation Ministry of Education and Science (grant MIREA 2.2.2.2.6546). 



[1] S. D. Drell and T.-M. Van, Phys. Rev. Lett. 25, 316 (1970). 

[2] S. Palestini et al, Phys. Rev. Lett. 55, 2649 (1985). 

[3] J. S. Conway et al, Phys. Rev. D39, 92 (1989). 

[4] E. L. Berger and S. J. Brodsky, Phys. Rev. Lett. 42, 940 (1979). 

[5] A. Brandenburg, S. J. Brodsky, V. V. Khoze, and D. Miiller, Phys. Rev. Lett. 73, 939 (1994). 

[6] A. Brandenburg, D. Miiller, and O. V. Teryaev, Phys. Rev. D53, 6180 (1996). 

[7] R. Bertini, O. Denisov, and S. Paul, Villars Meeting of the CERN SPSC, September 22-28, 

2004, Villars, Switzerland; Eol "Drell-Yan program at COMPASS", COMPASS Collaboration, 

in preparation. 

[8] A. P. Bakulev, S. V. Mikhailov, and N. G. Stefanis, Phys. Lett. B508, 279 (2001); B590, 
309(E) (2004). 

[9] V. L. Chernyak and A. R. Zhitnitsky, Phys. Rept. 112, 173 (1984). 
[10] A. V. Efremov and A. V. Radyushkin, Phys. Lett. B94, 245 (1980). 
[11] A. V. Efremov and A. V. Radyushkin, Theor. Math. Phys. 42, 97 (1980). 
[12] G. P. Lepage and S. J. Brodsky, Phys. Rev. D22, 2157 (1980). 



17 



[13] E. L. Berger, Z. Phys. C4, 289 (1980). 

[14] L. Del Debbio, Few Body Syst. 36, 77 (2005). 

[15] M. Gockeler et al, Nucl. Phys. Proc. Suppl. 161, 69 (2006), talk in the Workshop on Light- 
Cone QCD and Nonperturbative Hadron Physics 2005 (LC 2005), Cairns, Queensland, Aus- 
traha, 7-15 Jul 2005. 

[16] V. M. Braun et al, Phys. Rev. D74, 074501 (2006). 

[17] A. P. Bakulev, S. V. Mikhailov, and N. G. Stefanis, Phys. Rev. D67, 074012 (2003). 

[18] A. P. Bakulev, S. V. Mikhailov, and N. G. Stefanis, Phys. Lett. B578, 91 (2004). 

[19] A. P. Bakulev, S. V. Mikhailov, and N. G. Stefanis, Phys. Rev. D73, 056002 (2006). 

[20] J. Gronberg et al, Phys. Rev. D57, 33 (1998). 

[21] A. Schmedding and O. Yakovlev, Phys. Rev. D62, 116002 (2000). 

[22] A. P. Bakulev, S. V. Mikhailov, and N. G. Stefanis, Annalen Phys. 13, 629 (2004). 

[23] A. P. Bakulev and A. V. Pimikov, Acta Phys. Polon. B37, 3627 (2006). 

[24] K. Hagiwara et al, Phys. Rev. D66, 010001 (2002). 

[25] N. G. Stefanis, Nuovo Cim. A83, 205 (1984). 

[26] L. Del Debbio, Few Body Syst. 36, 77 (2005). 

[27] A. P. Bakulev, K. Passek-Kumericki, W. Schroers, and N. G. Stefanis, Phys. Rev. D70, 033014 

(2004); D70, 079906(E) (2004). 
[28] T. Gehrmann and W. J. Stirling, Z. Phys. C65, 461 (1995). 
[29] M. Gliick, E. Reya, and A. Vogt, Z. Phys. C67, 433 (1995). 
[30] C. S. Lam and W.-K. Tung, Phys. Rev. D21, 2712 (1980). 

[31] O. V. Teryaev, in Proceedings of the XI Advanced Research Workshop on High Energy Spin 
Physics (DUBNA-SPIN-05), September 27-October 1, 2005, Dubna, Russia, edited by A. V. 
Efremov and S. V. Goloskokov (JINR, Dubna, 2006), pp. 171-174. 

[32] B. Pire and J. P. Ralston, Phys. Rev. D28, 260 (1983). 

[33] R. D. Carlitz and R. S. Willey, Phys. Rev. D45, 2323 (1992). 

[34] H. Yokova. [arXiv:0705.24'8n [hep-ph]. 



18 



